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ABSTRACT 

The ~800-year old pulsar J1846— 0258 is a unique transition object between rotation-powered pulsars 
and magnetars: though behaving like a rotation-powered pulsar most of the time, in 2006 it exhibited 
a distinctly magnetar-like outburst accompanied by a large glitch. Here we present X-ray timing 
observations taken with the Rossi X-ray Timing Explorer over a 2.2 year period after the X-ray 
outburst and glitch had recovered. We observe that the braking index of the pulsar, previously 
measured to be n = 2.65 ± 0.01, is now n = 2.16 ± 0.13, a decrease of 18% ± 5%. We also note a 
persistent increase in the timing noise relative to the pre-outburst level. Despite the timing changes, 
a 2009 Chandra X-ray Observatory observation shows that the X-ray flux and spectrum of the pulsar 
and its wind nebula are consistent with the quiescent levels observed in 2000. 

Subject headings: pulsars: general — pulsars: individual (PSR J1846— 0258) — supernovae: individual 
(Kes 75) — X-rays: stars 
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1. INTRODUCTION 

PSR J1846-0258 is the very young (-800 year) X- 
ray pulsar at the center of the supernova remnant Kes 
75 (jGotthelf et all I2OOO0 . The pulsar has a spin pe- 
riod of 326-ms and usually exhibits properties common 
to those of rotation-powered pulsars (RPPs), including 
powering a bright pulsar wind nebula (PWN). This pul- 
sar is notable for having an unusually high spin-down- 
inferred magnetic field {B = 5 x 10^'^ G), and is one of 
the few with a measured braking index {n = vv jv^ = 
2.65 ± 0.01, where v is the spin-frequen cy, an d v and 
V are its derivatives; JLivingstonc ct al.j 120061 ). Mea- 
sured braking indices fall in the range 1.4 < n < 2.84 



(iLvne et al.lll993l[T996l; [Livingstone eralll2005allbl . l2Q0l 



iWeltevrede et al.ll2010D , all less than the canonical value 
of n = 3 pre dicted for mag netic dipole radiation in a 
vacuum (e.g.. lGunn fc Ostnk cr 1969). Explanations for 
rt < 3 include an increasing rnagnet ic moment (e.g., 
iBlandford fc Romanil [19881 lLvnell2004D or the effects of 
magnetospheric plasma on the spin-do wn torque (e.g., 
iHarding et "all 119991: Ifcameret al.|[2006l ). 

Against expectations, PSR J1846— 0258 exhibited dis- 
tinctly magnetar-like behavior in 2006 May- July, when 
it showed sev eral X-ray bursts, an X-ray flux increase 
(jGavriil et al.|[2008f ). a sizable r otational glitch with re- 
markable "overshoot" rec overy (jLivingstone et aLl 120101 : 
iKuiper fc HermsenI I2009D. as well as spectral changes 
(IKumar fc Safi-Harbl 12001 IGavriil etHlIM)! INg et al.l 
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l2008f ). This pulsar, evidently an RPP/magnetar transi- 
tion object, presents a unique opportunity of exploring 
the long-term relationship between magnetic activity and 
neutron star spin-down. 

Bursts of X-rays and variable X-ray flux, as observed 
in most magnetars and PSR J1846— 0258, are proposed 
to originate from small- or large-scale reorganiz ations of 
the magnetic held (e.g.. [Thompson et al.l [200^ . As the 
magnetic held is intimately connected with the temporal 
evolution of pulsars, comparing timing behavior before 
and after an episode of magnetic activity in a neutron 
star could provide important insight into the physics of 
neutron star magnetospheres. None of the magnetars 
have measured braking indices, despite their young ages 
(tc ~ 1 — 100kyr]3- This is because of a signiflcant level 
of low- frequency timing noise dCavriil fc Kaspi 2 004) and 
in several cases, large glitches (jDib et al.l 120081 and ref- 
erences therein), which prevent a measurement of n in 
every known magnetar. By contrast, PSR J1846— 0258 
rotates relatively steadily, allowing for a measurement of 
n with regular monitoring observations with the Rossi X- 
ray T iming Explorer (RX TE). Between its discovery in 
1999 fjGotthelf et al.ll2000f ) and the outburst and glitch 
in 2006, it experienced a level of timing noise typical 
of young RPPs along with one small glitch i^t^v jv ~ 
2.5 X 10~^ near MJD 52210) and one small candidate 
glitch {b.vjv < 5 X 10"® near MJD 52910). 

Also available from regular RXTE monitoring is the 

^ See1http://www. physics. mcgil l.ca/~pulsar/magnetar/main.html| 
for a catalog ot known magnetars. 
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pulsed flux o f PS R J1846-0258. As reported by 
iGavriil et alj ()2008[ ). the pulsed flux had returned to 
its quiescent value roughly two months after the 2006 
May outb urst and no further flux variations have been 
observed ([Livingstone et all I2010D . However, because 
RXTE is a non-focusing X-ray telescope, no information 
about the total flux or the phase-averaged spectrum are 
available from these data. To measure these quantities 
and to confirm the RXTE results, a focusing X-ray in- 
strument like the Chandra X-ray Observatory is required. 

Indeed, Chandra observations of the pulsar and PWN 
revealed fiux a nd spectral change s at the time of the 
2006 outburst (IGavriil et al.l 120081: iKumar fc Safi-Harbl 
120081 iNg et al.ll2008[ ). The pulsar's total flux rose con- 
siderably, but equally as interesting was its change in 
spectrum. The quiescent spectrum of PSR J1846— 0258 
is much like that of other young, high spin-down lu- 
minosity RPPs: a simple power law. Chandra obser- 
vations showed that while in outburst, the spectrum 
softened significantly such that it became reminiscent of 
those observed from Anomalous X-ray Pulsars (AXPs), 
namely, well described by a power law with an addi- 
tional thermal component. In addition, the superior 
angular resolution of Chandra allowed detailed obser- 
vations of the PWN, showing that a marginal increase 
in flux may have o ccurred between 2000 and immedi- 
ately post-outburst (jKumar fc Safi-Harbl [20081 ). The ef- 
fect of magnetar-like outbursts on PWNe is an open ques- 
tion, as none of the bona fide magnetars p ower nebulae 
(altho ugh see iVink fc Boiiibal ([2009i) and iTiengo et all 
([20101 ) for discussion of a possible PWN surr o undin g the 
AXP IE 1547.0-5408). IKumar fc Safi-Harbl (120081) sug- 
gested a causal relation between a possible increase in the 
PWN flux observed in 2006 and t his or p ast magnetar- 
like ou tbursts of the pulsar, while iKargal tsev & Pavloa 
(|20Q3) suggested that the PWN was over-luminous com- 
pared to those of other young pulsars, perhaps owing to 
previous, unseen magnetar o utbursts. However, a revised 
distance estimate of ~6-kpc ([Leahy fc Tian"2008'). rather 
than the previously claimed 19 kpc (Becker fc Hclfand 
I1984J ). reduces the implied nebular X-ray efficiency to 
77 = Lpwn/E ~ 0.02. While this is still large, it is sim- 
ilar to that observed from the Crab, hence need not be 
powered by previous outbursts. Thus any long-term ef- 
fect of the 2006 outburst on the PWN could help clarify 
this point. 

In this paper, we report on 2.2 year of RXTE timing 
observations of PSR J1846— 0258 in the post-magnetic 
activity, post-glitch recovery period. We perform phase- 
coherent and partially phase-coherent timing analyses af- 
ter the 2006 glitch had largely recovered, and report a 
post-burst measurement of n = 2.16 ± 0.13, smaller than 
the pre-outburst value at the 3.8cr level. We also quan- 
tify the increase of the timing noise observed over the 
bursting episode and discuss the implications of these 
observations. In addition, we report on 2009 Chandra 
observations of the pulsar and associated PWN which 
show that the pulsar and PWN fiux and spectra are con- 
sistent with their pre-outburst values. 

2. RXTE OBSERVATIONS AND ANALYSIS 

RXTE observations of PSR J1846-0258 were made us- 
ing the Proportional Counter Array (PCA: [Tahoda et al.l 
[r996; Ja hodact al. 2006,). The PC A consists of an ar- 



ray of five collimated xenon/methane multi-anode pro- 
portional counter units (PCUs) operating in the 2 - 
60 keV range, with a total effective area of approximately 
6500 cm2 and a field of view of ~ 1° FWHM. 

The RXTE data of PSR J1846-0258 are unevenly 
spaced over 11 years from 1999 April 18 through 2010 
Aprfl 22 (MJD 51286 - 55308). Observations taken be- 
tween 1999 Aprfl 18 and 21 (MJD 51286 - 51289) are 
excluded because they cannot be unambiguously phase 
connected with the rest of the data. Data from 2000 
January 31 to 2008 December 10 (MJD 51574 - 54810) 
were red uced and analyzed previo usly and details can be 
found in iLivingstone et al.l ([20060 and iLivingstone et al.l 
([2010D . Data taken between 2009 January 27 and 2010 
Aprfl 22 (MJD 54858 - 55308) are described here for the 
first time. 

Data were collected in "GoodXenon" mode, which 
records the arrival time (with l-fj,s resolution) and en- 
ergy (256 channel resolution) of every unrejected event. 
Typically, two to three PCUs were operational during an 
observation. We used the first Xenon layer of each op- 
erational PCU and extracted events in channels 4-48 
(approximately 2 - 20keV), as this produces good qual- 
ity profiles for individual observations. 

Observations were downloaded from the HEASARC 
archivcl^ and data from each active PCU were merged 
and binned at (1/1024) s resolution. Photon arrival times 
were converted to barycentric dynamical time (TDB) at 
the solar system barycenter using the J2000 source posi- 
tion RA^^18^W^2O4±0!01, deck = -02°58'30.1"± 
0.2" ([Helfand et al.ll2003D and the JPL DE200 solar sys- 
tem ephemeris. 

The phase-coherent ephemeris from ILivingstone et aLl 
([2O IOI) was used to fold the time series for each new obser- 
vation with 16 phase bins. Resulting profiles were cross- 
correlated with a high significance template profile, pro- 
ducing a single Time-Of-Arrival (TOA) for each obser- 
vation. This process produced 265 TOAs between MJD 
51574 and 55308 with a typical uncertainty of ~ 10 ms 
(~3% of the pulse period). The TOAs were fitted to 
a timing model using the pulsar timing software pack- 
age TE MPCEl. Further details of the timing analysis are 
given in ILivingstone et al.l ([2006[) . 

3. TIMING ANALYSIS AND RESULTS 

3.1. Phase- coherent Timing Analysis 

In order to make a significant measurement of a deter- 
ministic value of iy and thus n, we restricted our phase- 
coherent timing analysis to MJD 54492 - 55308 (2008 
January 27 ~ 2010 April 22), because earlier observations 
are highly contamin ated by glitch recovery a nd timing 
noise, as discussed in ILivingstone et al.l ([2010f ). 

We obtained a single phase-coherent timing solution 
(with no phase ambiguities) fitting only v and iy for the 
100 TOAs during this time period, shown in the top 
panel of Figure [TJ The residuals show a very signifi- 
cant phase contribution from a second frequency deriva- 
tive (i.e., ~15 phase turns over roughly 2 years). We 
therefore added i) to our phase-coherent fit, and the re- 
sulting residuals are shown in the middle panel of the 
figure. The fitted i> corresponds to a braking index of 

* http://heasarc.gsfc.nasa.gov/docs/archive.html 
^ http://www.atnf.csiro.au/research/puIsar/tempo/ 
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n = 1.888 ± 0.002. As visible in the middle panel of 
Figure [U significant timing noise remains in the data; 
the timing residuals are not Gaussian distributed. As a 
result, the formal la uncertainty on n from this global 
phase-coherent fit underestimates the true uncertainty. 
Spin parameters from this fit are given in Table [TJ 

We fitted higher order frequency derivatives to 
"whiten" the phase residuals, a common procedure to 
lessen the conta minating effect of timing noise on fitted 
parameters (e.g.. lKaspi et al.iri994l ). Fitting 12 frequency 
derivatives (the maximum possible given the current ma- 
chine precision) removes the majority of the timing noise, 
though does not render the phase residuals entirely Gaus- 
sian, as shown in the bottom panel of Figure [1] Table [2] 
shows the variation of n as derivatives are fitted, as well 
as x^i^- The value of n varies between 1.89 and 2.95 as 
higher order derivatives are fitted, without converging to 
a single value, rendering the true value of n ambiguous 
from this analysis. Nevertheless, the range of measured n 
values from this analysis is relatively narrow: the timing 
noise has increased and clearly contaminates, but does 
not completely dominate v. In cases where a parame- 
ter is dominated by a noise process, it can be of several 
order s of magnitude lar ger and often of the wrong sign 
(e.g., iHobbs et al.ll201Cil ). neither of which are seen here. 
In cases where timing noise contaminates a measurement 
of a deterministic parameter but may not dominate, as 
for these PSR J1846— 0258 data, using a partially coher- 
ent timing analyis can be useful to find the true value. 

3.2. Partially Coherent Timing Analysis 

To mitigate the effects of timing noise, we performed a 
partially coherent timing analysis. We created 48 short 
phase-coherent timing solutions spanning from 2000 to 
2010 (MJDs 51574-55308). For each short timing solu- 
tion, we fit only v and i). The time span included in each 
subset was determined from the requirement that the re- 
duced x^ of the fit was ~1, and that no red noise- like 
structure was visible in the data. In addition, we allowed 
the data included in each timing solution to overlap by 
^1/2, where sufficiently dense sampling was available. 
This improves coverage, which is of particular impor- 
tance while the timing noise level is very high. Figure [5] 
shows the resulting v measurements spanning from 2000 
to 2010. This analysis excludes observations taken be- 
tween 2006 May 31 and 2007 January 27, when glitch re- 
covery and timing noi se prevented a coh e rent t iming so- 
lution. As discussed in [Livingstone et al.l (|2006() and con- 
firmed by our current analysis, from 2000 - 2006 May, i> 
increased very regularly, except at a small glitch in 2001. 
The steady increase in v corresponds to a braking index 
of n = 2.65 ± 0.01, as measured from a phase-coherent 
analysis of these data. The large glitch (visible as a sud- 
den decrease in z> by ~3% in Figure [5]) and the increase 
in timing noise had largely recovered by the beginning of 
2008, as shown in the figure. 

To obtain a post-burst measurement of n, we ignored 
ah timing data prior to MJD 54492 (2008 January 27), 
where the aforementioned glitch recovery and timing 
noise dominate. We performed a weighted least-squares 
fit to 16 i> measurements spanning MJDs 54492 - 55308 
(2008 January 27 - 2010 April 22, shown in the inset 
of Figure 2). Given the large scatter in the post-burst 
z> measurements, and the known effects of timing noise 



on these data, it is likely that the formal uncertain- 
ties significantly underestimate the true uncertainties. 
Thus, to better estimate the uncertainty on i/, we em- 
ployed a bootstrap error analysis, helpful when formal 
uncertainties may underestimate the true uncertainties 
(Efron 1979i). and previously us ed for this same pur- 
pose in [Livingstone et all ()2005al lb[) . This results in i> = 
3.13(19) X 10-21 s-3, corresponding to n = 2.16 ± 0.13, 
where the uncertainty from the bootstrap estimate is 
larger than the formal uncertainty by a factor of ~2.4. 
This new measurement of n is smaller than the pre- 
outburst value of n = 2.65±0.01 at the 3.8ct level (or 9.1cr 
if only the formal uncertainties are considered). Thus, 
the braking index decreased by An = —0.49 ± 0.13, fol- 
lowing the period of magnetar-like activity in 2006. This 
is the first observed significant measurement of a change 
of a braking index. 

To further confirm a change in n, we performed an 
identical partially coherent timing analysis on the pre- 
outburst measurements of v. Because of the small glitch 
and candidate glitch, the data were separated into three 
large subsets. The first subset (prior to the small glitch 
near MJD 52210) contains only three measurements of 
z>, resulting in n = 2.63 ± 0.04 (where uncertainties for 
this measurement are the formal uncertainties because 
a bootstrap error analysis cannot be performed with no 
additional degrees of freedom). The second subset lies 
between the small glitch and a 78 day gap in the data 
(containing a candidate glitch). In this data subset, we 
performed a weighted least squares fit and a bootstrap 
error analysis on seven measurements of v. This re- 
sulted in ri = 2.61 ± 0.07. We repeated this analysis 
for the eight measurements of u prior to the magnetar- 
like outburst, resulting in n = 2.68 ± 0.03. The three 
values of n are in agreement with each other and with 
the value obtained from a fully phase-coherent t iming 
analysis (n = 2.65 ± 0.01; iLivingstone et al.l 120061 ). All 
pre-out burst measurements of n are systematically larger 
than that measured post-outburst. Note that the un- 
certainties on the two measurements of n from partially 
coherent analyses (with bootstrap uncertainties) to pre- 
outburst data are smaller than the uncertainty for the 
post-outburst value of n, despite similar data spans fit- 
ted in each case. This is indicative of an increase in the 
timing noise post-outburst, discussed further below. 

A complicating factor in timing some magnetically 
active neutron stars is that pulse profile changes of- 
ten accompa.ny ra diative changes and/or glitches (e.g., 
iKaspi et al.l 120031 ). This can affect measured timing 
parameters and must therefore be quantified. It has 
previously been reported that no significant changes 
in the pulse profile were detected during the ou tburst 
(jLivingstone et al.l 120101: iKuiper fc HermsenI [20091 ) . We 
further verified that the pulse profile in the ~2 year pe- 
riod immediately prior to the outburst is not statistically 
different from the summed profile from the ^^2 year of 
data used to measure n, shown in Figure |3l 

In order to further analyze the timing noise contami- 
nating these data, we performed a second partially co- 
herent timing analysis, this time with each data subset 
having z/, i/, and iy fitted. The same conditions oix^u ~ 1 
and Gaussian-distributed residuals, were applied to de- 
termine the length of each subset, and once again, the 
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period between 2006 May 31 and 2007 January 27 was 
excluded owing to the lack of a coherent timing solution. 
Six values of i) were obtained before the outburst, while 
nine values of i) were obtained after the outburst, shown 
in Figure |4l The measurements of v before the outburst 
indicate the regular rotation of the pulsar during this pe- 
riod, while the single value of i) above the average (visi- 
ble in the inset) occurs directly after the candidate glitch 
near MJD 52910, providing the best evidence for a glitch 
during this period. As visible in the figure, the value of v 
changed dramatically immediately after the outburst, to 
a maximum of ^200 times the quiescent value, as well as 
varying in sign, indicating a dramatic increase in timing 
noise during the period of glitch recovery. The magnitude 
of V decays as the glitch recovers during 2007. The inset 
of Figure |4] shows the variation of i^ on a smaller scale. 
While the effects of glitch recovery and timing noise have 
subsided substantially by 2008, the post-outburst varia- 
tion of i) remains larger than its pre-outburst behavior. 

3.3. Timing Noise 

Qualitatively, the timing noise in the 2.2 year period 
used to obtain the post-burst measurement of n is larger 
than that observed prior to the outburst, though much 
smaller than in the initial aftermath of the outburst, 
when no phase- coherent timing solution was possible. 
One measure of the change in timing noise can be found 
by comparing the RMS residuals from MJDs 54492 - 
55308 with those in a similar time span before the out- 
burst. Fitting V, !>, and i) phase coherently to the 2.2 
year segment of data spanning MJD 53086 - 53879, just 
before the outburst, results in a timing solution with 
RMS residuals of 11.4-ms (0.035 periods), a factor of 
--5.5 smaller than the RMS residuals from MJDs 54492 - 
55308, of 63.6ms (0.19 periods). This shows that the 
timing noise post-burst is significantly larger than be- 
fore the magnetar-like outburst. This is refiected in the 
new measurement of n by the increased uncertainty com- 
pared to that on the pre-outburst value of n for similar 
measurement baselines. 

Another measure of the increase in timing noise is the 
time span that can be included in each partially coherent 
measurement of v as shown in Figure [5] When the pulsar 
is less noisy, more data can be included in each short tim- 
ing solution while satisfying the conditions that x^,^ ^ 1 
and that no red noise-like structure remain in the data. 
The pulse profile and pulsed flux are steady, important 
because variability in either could cause changes in TOA 
uncertainties, and thus affect the time span for each v 
measurement. From 2000 until 2006 May, on average, 
each measurement of v was obtained with data spanning 
111 ± 26 days, while in 2007 (when glitch recovery was 
still a significant effect), each measurement of v spanned 
an average of 33 ± 20 days. From 2008 to 2010, the av- 
erage time span for each measurement was 68 ± 16 days. 
Thus, nearly four years post-outburst, the pulsar remains 
noisier than prior to the outburst. 

A well known measure of timing noise is the Ag param- 
eter, defined as the contribution to the rotational phase 
of the pulsar from a measurement of D over a period of 
10^ s assuming that i) is e ntirely dominated by timing 
noise ([Arzoumanian et al.l 11994). This parameter is of 
limited value for a pulsar where i) is dominated instead 
by secular spin-down, where most of the phase contribu- 



tion from i) is due to magnetic braking or another deter- 
ministic spin-down mechanism. To quantify the change 
in timing noise observed in PSR J1846— 0258, we define 
a similar parameter which quantifies the contribution to 
the rotational phase from the measurement of the third 
frequency derivative, v , over a time span of ^ 2.5 x 10^ s. 
The time span is optimized for this particular pulsar, as 
it is the approximate amount of time required to obtain 
a significant measurement of v , while allowing several 
measurements to be made given the available data span. 
Thus, in analogy with the Ag parameter we define: 
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We measured the A- parameter for each ^ 2.5 x 10^ s 
segment of data where a phase-coherent timing solution 
was available, and have shown the results in Figure [5j 
The value of A- increases dramatically after the 2006 
outburst, after which it decays, but by 2010 has not re- 
turned to the pre-outburst quiescent level. 

4. CHANDRA OBSERVATIONS AND ANALYSIS 

PSR J1846— 0258 was previously observed with the 
Chandra X-ray Observatory in 2000 and 2006 with 37-ks 
and 155-ks exposures, respectively. Serendipitously, the 
latter observation was carried out seven days after the 
outburst. We obtained a new Chandra imaging obser- 
vation with a 44.6-ks ACIS-S exposure (ObsID 10938) 
on MJD 55053 (2009 August 10), over three years after 
the outburst. The observation was taken in 1/8 subar- 
ray mode, which gives a short frame time of 0.4-s, re- 
sulting in negligible (< 3%) pileup of the pulsar counts. 
All the data reduction and analysis were performed with 
CIAO 4.10. Figure [B] shows the exposure-corrected im- 
ages of PSR J1846-0258 and the associated PWN in 
the 1 - 7keV energy range. While the PWN exhibits 
no obvious change in the overall morphology between 
the three observations, time variabilities are observed in 
some small-scale features. The most prominent area of 
variability is the northern clump located ^ 7" northeast 
of the pulsar. We extracted the count profiles of the 
clump from the exposure-corrected images using 6" wide 
boxes (as indicated in Figure |6]), and show the results for 
each epoch in Figure [T] The clump morphology evolved 
from a single peak in 2000 to a double peak in 2006, and 
back to a fainter single peak in 2009. However, there is 
no evidence for bulk motion of the clump, as is visible in 
the figure. 

For the spectral analysis, we extracted the spectrum 
of the entire PWN excluding the central 3" radius, and 
fitted it with an absorbed power-law model. The col- 
umn density was fixed at A^'h = 4.0 x lO^^cm"^ dur- 
ing the fit, to provide a di r ect comparison with prev i- 
ous studies (|Ng et all 120081: iKumar fc Safi-Harbl [2008[) . 
In the 2009 observation, we found a photon index of 
r = 1.90 ± 0.03 and an absorbed flux of fti-w = 
(1.37 ± 0.05) X 10""ergs-icm-2 in the 0.5 - lOkeV 
range (hereafter, quoted uncertainties are 90% confidence 
levels). This is fully consistent with the PWN spec- 
trum in 2000, which has F = 1.88 ± 0.03 an d /(f.^'_io = 
(1.39 ±0.02) X 10-"crgs-icm-2 (INg et al.l [20081 . 

^ |http : //cxc ■ harvard ■ edu/clao4 .i./\ 
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For the pulsar itself, counts were extracted from a 
2"radius aperture and grouped with 70 counts in each 
spectral bins. We followed Ng et al. (2008) to account for 
the nebular contamination using a power-law model with 
fixed r = 1.9 and fi^f_io = 1.0 x lO-^^gj-gs"! cm^^. 
As in previous studies, we started with a power-law plus 
blackbody model, but found that the latter component is 
only marginally detected, statistically insignificant. The 
pulsar spectrum in 2009 is adequately fitted by an ab- 
sorbed power-law model with F — 1.1 ±0.1, and we derive 
an upper limit on the blackbody temperature of 0.25 keV 
at the 90% confidence level, significantly lower than the 
value 0.9±0.2keV measured during the 2006 outburst. 
Table[3l compares the best-fit spectral parameters among 
the three epochs, indicating that the pulsar fiux in 2009 
is at a similar level as in 2000, much lower than that in 
2006. For completeness, we also fitted a two blackbody 
model to the 2009 Chandra observation. While the good- 
ness of fit is statistically similar to that obtained when fit- 
ting a power-law model, the two-blackbody model gives 
an unphysically high temperature of 1.6± 0.1 keV for the 
hotter thermal component. In addition, there is evidence 
for a continuous pow er-law component from 1 to 300 keV 
(|Kuit)er fc Hermsen 2009), rendering the two blackbody 
model of limited interest for this pulsar. 

5. DISCUSSION 
5.1. Timing Noise and the Braking Index 

We have observed a decrease in the braking index of 
PSR J1846-0258 from n = 2.65±0.01 to n = 2.16±0.13, 
corresponding to An — —0.49 ± 0.13, or a decrease of 
18% ± 5%. The change in n was accompanied by an 
increase in the timing noise of the pulsar, which remains 
larger than the pre-outburst level, nearly four years after 
the glitch and outburst on 2006 May 31. Previous long- 
term observations of n in young pulsars have shown that 
they are remarkably steady. In PSR B1509— 58, timing 
observat ions over 21 years show that n varies by only 
~ 1.5% (jLivingstone et al] |2005b'). whi le the Crab pulsa r 
exhibits variations on the order of 5% (jLvne et al.lll993[ ). 

There are two possible ways to interpret the measure- 
ment of An = -0.49 ± 0.13. The first is that the true 
n decreased by a significant amount after the 2006 out- 
burst. The second is that the increased timing noise is 
causing an apparent decrease in n. We discuss each of 
these interpretations next. 

5.1.1. Variable Braking Index 

If the true value of the braking index changed perma- 
nently at the time of the magnetar-like outburst, what 
could be the physical cause of such an effect? From the 
spin-down law derived assuming magnetic dipole braking 
(where / is the moment of inertia, R is the neutron star 
radius, an d a is the angle be tween the spin and magnetic 
axes; e.g.. lGunn fc Ostriked[l969.' ). 



~3 J? ' 



(2) 



we can infer that d^I/dt^ > CQ, da/dt > 0, or dB/dt > 
will cause n < 3. 

'^ Note that the different relationship between I and An than for 
B and a arises from correctly including the variability of / in the 



While it is hard to imagine a physical situation caus- 
ing an accelerated growth or decay in /, varying values 
of a or B have been considered in the past. Counter- 
alignment of the magnetic field (i.e., an increasing a) 
results in n < 3, even if An = 0. A sudden increase in 
the rate of change of a would produce An < 0. How- 
ever, this is difficult to invoke for the observations of 
PSR J1846— 0258 because the pulse profile shows no vari- 
ation over the relevant time period (e.g.. Figure [3]). A 
small change in a could still be possible if our line of 
sight is crossing the center of the pulsar beam, as a large 
change in the pulse profile may not be required from a 
small change in field orientation. However, this is hard 
to reconcile with the lack of detected radio pulsations 
from the source (| Archibald et al.l 120081) . which is typi- 
cally interpreted as our line of sight missing the magnetic 
pole, from where the radio emission is thought to origi- 
nate, and crossing only the wider X-ray beam. Alterna- 
tively, radio emission from PSR J1846— 0258 may be sup- 
pressed as a result of its larg e magnetic field, as sugge sted 
for most known magnetars ([Baring fc Hardin^ll998[ ) , al- 
though radio emission has been detected from several 
pulsars with higher inferred magnetic dipole strengths 
than observed from PSR J1846-0258. Nevertheless, if 
radio emission were suppressed in PSR J1846— 0258, then 
our line of sight may cross the center of the X-ray beam, 
allowing for a small change in a without accompanying 
profile changes. This point cannot be resolved, however, 
without further details of the geometry of the system. 

An increase in the magnitude of B, without a change 
in the orie ntation of the field could also cause a de- 
creas e in n (jBlandford fc Romanilll988l : iBlandford et al.1 
119831 ). Making all the assumptions of perfect dipole 
spin-down but allowing dB/dt > 0, a braking index of 
n = 2.65 implies a time scale of growth for the magnetic 
field of ^8000-year, while n == 2.16 shortens the growth 
timescale to ~3500-year. The possibility of magnetic 
field growth is intriguing given the magnetic activity ob- 
served prior to the change in n. iBhattacharva fc Sonil 
IJOO?) suggest that magnetar-strength fields emerge over 
a period of time, as shielding currents dissipate. If the 
effective B is currently in such a period of growth, the 
smaller value of n could result (jLynd |200^ . However, 
this picture lacks a description of why the B-growth 
timescale should have increased so rapidly. 

The standard spin-down law (Equation [2|) is a ma- 
jor idealization even for rotation-powered pulsars; for 
magnetars, the picture is almost certainly much more 
complicated. According to one version of the magnetar 
model, there is a one-to-one relationship between n and 
the large-scale twist angle between the north and south 
hemi spheres of magnetic field, A0n-s (IThompson et al.) 
I2002f ). Here, a decrease in n would imply an increase 
in the twist angle. A braking index of n = 2.65 corre- 
sponds to a twist angle of A0 ~ Irad, while n ~ 2.16 
implies a twist angle of A^ ~ 2rad. If the magnetic 
field remains approximately constant before and after the 
outburst, such an increase in the twist angle should be 



spin-down luminosity, E, when deriving the above spin-down law. 
The full form of the spin-down luminosity is -B = AttIuu + 4-ku^I, 
where only the first term is considered if / is constant. Using the 
full form of E in the derivation of the spin-down law leads to the 
dependence of cPl /dt^ on An. 
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accompanied by an increase in the X-ray luminosity of 
'^ 50%, whereas no persistent increase in Lx is observed 
in PSR J1846— 0258 as evidenced by the consistent flux 
value measured with Chandra in 2000 and 2009. Fur- 
thermore, while the above mod el assu mes a global mag- 
netic field twist, Beloborodovl ()2009( ) argues that such 
global self-similar twists do not present a viable expla- 
nation for magnetar behavior. Instead, he suggests that 
magnetospheric currents are confined to narrow regions 
on the most extended field lines, which are not respon- 
sible for the bulk of magnetar X-ray emission. If this 
picture is correct, the spin-down of the star can vary 
without accompanying changes in X-ray luminosity. In 
addition, this model predicts that an increase in spin- 
down torque (though not necessarily monotonic) should 
occur sometime after a radiative event, and should even- 
tually return to the pre-outburst torque value. Qualita- 
tively, this seems to present an explanation of the ob- 
served timing variability in PSR J1846— 0258; however, 
it provides no quantitative prediction for n. A variable n 
provides an u nambiguous tes t of the theory of spin-down 
presented bv iMelatosI (|1997t ). He posits that the radius 
relevant to neutron star spin-down is not the point-like 
neutron star radius, but a somewhat larger "vacuum ra- 
dius" where field aligned fiow breaks down. This radius is 
large enough that the system can no longer be treated as 
a point dipole, resulting in modifications to the standard 
spin-down predictions. In the context of this model, a 
measurement oi v, v, and a uniquely predicts n. While 
there is no estimate of a for PSR J1846— 0258, there 
are now two measurements of n so that the theory can 
be tested. Given the observed change in n, and assum- 
ing that a is stable over the magnetar- like outburst, this 
theory predicts that the magnetic field should have in- 
creased by a factor of ^6. This is not observed, however, 
as the magnetic field estimate has increased by just 0.3% 
compared with the pre-outburst value, in contradiction 
to this theory. Alternatively, if B is roughly constant 
as observed, and inst ead a we re allowed to change at 
the time of the event, iMelatod (1997) predicts that the 
angle between the spin and magnetic axes would have 
changed from ~9° to ~4°. Such a change might be vis- 
ible as differences in the pulse profile for some lines of 
sight, but cannot be excluded given the lack of observed 
profile changes. 

Another possible explanation put forth for a static 
measurement of n < 3 is that magnetic field 
lines ar e deformed due to plasma in the magneto- 
sphere (jBlandford fc Romanil 119881 ). A sudden in- 
crease in the amount of plasma in the magnetosphere 
of PSR J1846— 0258, perhaps injected at the time of 
the magnetar-like outburst, could cause n to decrease. 
The best evidence for a plasma-filled magnetosphere af- 
fecting pulsar spin-down is found from the "intermit- 
tent" radio pulsar, PSR B1931-I-24, which has dramatic, 
quasi-periodic changes in i> corre lated with radio pu l- 
sations which turn on and off (jKramer et al.l I2006D . 
iHarding et al.l (1999) propose that spin down can arise 
from a combination of magnetic dipole radiation and 
wind losses. An increase in losses from the wind relative 
to dipole radiation will manifest a s a smaller value of n . 
The spin-down formula given by IHarding et al.l (|1999( ) 



implies a braking index of 



2vL„^/'^BR^ 
n = 6 , 



(3) 



where Lp is the kinematic luminosity of the wind, which 
can, in turn, be estimated with a measurement of n. For 
PSR J1846-0258, a change in n from 2.65 to 2.16 corre- 
sponds to nearly an order of magnitude increase in the 
persistent particle luminosity from Lp ~ 1 x lO^^ergs"^ 
to ip ~ 6 X lO'^^ergs^^. Because the outfiowing parti- 
cles travel near the speed of light, the additional particles 
would populate the PWN and be observable as a factor 
of ~6 increase in the PWN flux in the 2009 Chandra ob- 
servation as compared to the observation in 2000. How- 
ever, no such fiux increase is detected (see Section 4). 
The lack of a PWN flux increase seems to refute the idea 
that an increase in wind losses is responsible for An < 0, 
however, a more rigorous derivation of the relationship 
between wind losses and spin down may provide further 
insight into this issue. 

Interestingly, because a change in the plasma condi- 
tions in the magnetosphere might also affect magneto- 
spheric torques, t his could poss ibly explain a n increase in 
timing noise fe.g.. lChen^ll987[ ). As shown bv lLvne et al.l 
(|2010f) . timing noise in some pulsars can he traced to 
magnetospheric fluctuations. For these radio pulsars, 
there is a correlation between torque variations and pulse 
shape changes. This is difficult to apply in the case of 
PSR J18 46— 0258 because no ra dio pulsations have been 
detected (jArchibald et al.l[2008[ ) , and there is no evidence 
fo r profile variability in the X-ray b and (see Figure [3]). 

iContopoulos fc Spitkovskvl (|2006D note that if the co- 
rotation radius of the magnetosphere, rco-rot , is less than 
the light cylinder radius, i?LC, owing to imperfect re- 
connection within the magnetosphere, 1 < n < 3 will 
result. If the co-rotation radius decreased at the time 
of the magnetar-like outburst from magnetic field lines 
opening, n would also decrease. They parameterize the 
relationship between Tco-rot and n as 



^co-rot 



^LC I — 



(4) 



If the initial spin-period of PSR J1846-0258 was Pq = 
Ims, the implied co-rotation radius prior to outburst 
was Tco-rot = 0.36i?LC when n = 2.65 ± 0.01, and 
would have decreased to 0.09i?LC post-outburst when 
n = 2. 16 ±0.13. Changes to the extent of the co-rotation 
region of the magnetosphere might be visible as changes 
in the pulse profile; however, no such changes have been 
observed from PSR J1846-0258. In the context of this 
model, however, it is impossible to rule out changes to 
the extent of the magnetosphere as no independent esti- 
mate of the initial spin period is available. In the context 
of this model, taking Pq = 1 ms as the lower limit on the 
birth spin period for PSR J1846— 0258, provides a lower 
limit on the co-rotation radius for any measured value of 
n. Furthermore, Pq = 1 ms provides an upper limit on 
the change in co-rotation radius at the two epochs, since 
this quantity decreases with slower birth spin periods. 

5.1.2. Timing Noise 
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The second possible interpretation of An < is that 
the true n is constant but remains masked by the high 
timing noise. The increase in timing noise could arise as a 
result of changes to the superfluid interior brought on by 
the unusual 2006 glitch or changes in the magnetosphere 
after the outburst. Though a bootstrap error analysis 
was employed in order to better account for the effect 
of the increase in timing noise, a definitive test is not 
possible. Continued timing observations may be able to 
solve this issue, if the timing noise level continues to de- 
crease and the new braking index remains steady. How- 
ever, it is also possible that the increased level of timing 
noise and the decreased braking index are connected, for 
example, via an increase in the magnetospheric plasma 
density. In that case, if the pre-outburst conditions are 
eventually reobtained, both the timing noise and braking 
index should relax to their pre-outburst values, rendering 
a temporary value of n < 2.65 ambiguous in nature. 

Fluctuation in pinned superfluid in the pulsar inte- 
rior is one of the possible causes of timing noise (e.g., 
lAlpar et al.l 119861 ). Thus, one possible explanation for 
the increase in timing noise is that significant changes 
were imparted to the neutron star interior at the time of 
the 2006 glitch and outburst. The glitch was followed by 
an unusual over-shoot recovery (Q ~ 8.7 on a timescale 
of Td ~ 127 days) and a permanent increase in the mag- 
nitude of ;> (w ith fractional magnitude Ai//i' ~ 0.0041; 
[Livingstone et al. 2010). This recovery is thus far unique 
among pulsar glitches and the origin and long-term con- 
sequences of such behavior are not well understood. By 
contrast, the permanent change in i) following the glitch 
is not unusual when compared to those measured after 
other glitches. In addition, we note that the change in ;> 
is not responsible for the observed decrease in n. Because 
the fractional increase in i) is three orders of magnitude 
smaller than the fractional change in i), it is the latter 
effect that dominates the change in n. 

In addition to the detected glitch recovery and perma- 
nent increase in ;>, non-monotonic variations in ;> were 
observed in the aftermath of the glitch. While glitch 
recovery (i.e., a temporary increase in i/) and discrete 
jumps in v accompanying a glitch are both established 
phenomena, to our knowledge, no other RPP has expe- 
rienced changes in timing noise similar to those observed 
from PSR J1846— 0258. However, variable spin-down 
torque has been observed in several magnetars. The 
AXP IE 1048.1—5937 twice experienced approximately 
year-long periods of rapid i> variations, i.e., sudden, but 
temporary increases in timing noise (jGavriil fc Kaspj 
120041 : iDib et ahl I2009D . Similar variations in v were ob- 
served in the transien t AXP XT E 1810-197 after its 
2003 outburst (Camilo et al.ll2007D . Thus, the observed 
change in timing noise in PSR J1846— 0258 can be inter- 
preted as yet another example of magnetar-like proper- 
ties from this RPP, even if the phenomenon is currently 
unexplained. 

5.2. Phase-averaged pulsar flux and the pulsar wind 

nebula 

The Chandra results show that the phase-averaged pul- 
sar fiux and spectrum in 2009 have returned to the qui- 
escent values observed in 2000. This agrees with the 
pulsed flux history from RXTE, where the initial fiux 
increase in 2006 was observed to decay exponentially 



with 1/e ~ 55 days (jGavriil et al.ll2008D. and has sinc e 
remained at the quiescent level (jLivingstone et al.|[2"010|) . 
Given that the 2009 Chandra data were taken over 1000 
days after the 2006 outburst, our findings are not unex- 
pected. The time variability of the small-scale features in 
the PWN could be attributed to magnetohydrodynamic 
instabilities in the fiow, similar to what has been ob- 
served in th e PWNe powered by the Vela pulsar and PSR 
B1509-58 (jPavlov et al.|[200l iDeLanev et al.ll2Q06l ) and 
is unlikely to be related to a sudden deposition of par- 
ticles at the time of the outburst, ^avriil et al. (200^ 
reported an energy release of ~ 5 x 10'*^((i/6kpc)^ erg 
(2 - 60keV) in the 2006 outburst. In comparison, the 
total energy released in the giant fiare fro m the magne- 
tar S GR 1806-20 in 2004 is 2 x 10^^ erg (jPalmer et al.l 
l2005fl. of which 4 x 10^'^ erg went into particle energy 
(jGaensler et al.ll2005| ). The resulting ratio between the 
particle energy and electromagnetic radiation is 2 x 10~^. 
If the same ratio holds for PSR J1846-0258, then the en- 
ergy of the injected particles would be only ^ 10'^^ erg. 
With a B-field strength of 15^ G in the PWN (Djannati- 
Atai et al. 2007), the synchrotron cooling timescale of 
a particle emitting at 5-keV is ~300-year. These parti- 
cles would have induced an X-ray fiux enhancement in 
the PWN of 2 X 10~^^ergs~^ cm~^. This is five orders 
of magnitude lower than the X-ray flux of the northern 
clump, and six orders lower than that of the entire PWN, 
too small to be detected. Moreover, if the particles travel 
isotropically with a typical post-shock fiow speed of c/3, 
the fiux would have already spread over the entire PWN 
of 10" radius. Therefore, we do not expect the 2006 out- 
burst to have any observable effects on the PWN. 

6. CONCLUSIONS 

The observed change in n after the magnetar-like out- 
burst in PSR J1846— 0258, if shown to be steady via 
ongoing timing observations, has important implications 
for the physics of neutron star spin-down. A decrease in 
n could have several origins, and a detailed theoretical 
framework is necessary for interpreting this observation. 

Most theoretical descriptions of a changing n require 
an accompanying persistent change in radiative behav- 
ior of the pulsar, while we observe neither pulse profile 
variability or persistent fiux enhancement. An increase in 
particle wind losses relative to dipole losses does not pro- 
vide a good description of An < for PSR J1846-0258 
because o f the lack of a persis tent increase in PWN lu- 
minosity ([Harding et al.iri999[ ). However, variability in 
magnetospheric plasma remains a promising avenue for 
future consideration, especially considering the recent 
report of variable spin-down rate correlated with radio 
pulse shape changes for several pulsars, which confirms 
a link b etween torque and emission properties in several 
pulsars (jLvne et al.l 120101 ). No variability in the X-ray 
pulse profile is detected in PSR J1846— 0258; however, 
small shape changes may be present but not observable 
in current data. In addition, short time-scale variability 
would not be detectable in RXTE observations which are 
typically from 1.5 to 2hr long. 

The timing noise in PSR J1846— 0258 is observed to 
be of a higher level than prior to the outburst. That 
is, even four years after the glitch and magnetic activ- 
ity, the pulsar is rotating less regularly than in its pre- 
outburst quiescent state. It is interesting to note, how- 
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ever, that the current timing noise in PSR J1846— 0258 
would not be unusual if observed in any young pulsar. 
Rather, it is the sudden change in the level of timing 
noise in PSR J1846— 0258 that is noteworthy. Since the 
timing noise is simply of a higher level and not other- 
wise different from that observed in other pulsars, the 
phenomenon cannot be used as a diagnostic of previous 
unseen magnetic activity in other pulsars. 

The observed decrease in n and increase in timing noise 
reported here may or may not be permanent. Regu- 
lar monitoring observations beyond the RXTE era may 
help to answer this question, as well as to search for 
future magnet ar- like X-ray outbursts and glitches from 
PSR J1846-0258. 

The 2009 Chandra observations of PSR J1846-0258 
show that the total flux and spectrum are consistent with 
the quiescent values observed in 2000. No significant 
variability was detected in the PWN that can be asso- 
ciated with the 2006 outburst, nor, given the energetics 
of the outburst, would any be expected. The variability 
of the PWN observed in the 2006 Chandra observation 
is most likely unrelated to the outburst, and is probably 
similar in origin to the variation in small-scale features 
seen in other PWNe. 
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Figure 1. Timing residuals of PSR J1846-0258 spanning MJDs 54492 - 55308 (2008 January 27 - 2010 April 22). The top panel shows 
residuals with u and u fitted. The middle panel shows residuals with u, 0, and i> fitted, while the bottom panel shows residuals with 12 
frequency derivatives fitted. 
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Figure 2. Evolution of the frequency derivative of PSR J1846— 0258 over ~10-years of RXTE timing observations. Measurements of u 
overlap by ~l/2 where sufficient data are available. The effect of the two confirmed glitches is visible in the figure, the first near MJD 52210 
as a small discrete jump in u and the second as the large decrease in u near MJD 53886 (note that the full effect of the glitch on i> is not 
shown here in the interest of making visible the smaller changes in u at other epochs, see Livingstone et al. (2010) for additional details 
about the large glitch). The inset shows measurements of u from MJDs 54492 - 55308, with the best-fit slope shown as a solid line and the 
ztlfT uncertainties (from the bootstrap analysis which better accounts for the increase in timing noise) shown as dotted lines. The slope 
corresponding to the pre-outburst n is shown as a dashed line (colored blue in the online version) , where uncertainties are roughly an order 
of magnitude smaller than those from the post-outburst era, so are not visible in the figure. 



• I — I 

m 
0) 



T3 
m 

3 
ex 

N 

• i-H 

a 

o 



Post-Outburst Observations of PSR J1846-0258 

0.5 1 1.5 



11 



T I I I I I I I I I I I I I I I I r 



Tl 








MJD 53554 

Exposure 304 ks_ 



MJD 54827 

Exposure 409 ks 






J I I L 



J I I L 



J I I L 



J I I L 







0.5 1 

Pulse Phase 



1.5 



Figure 3. Pulse profile of PSR J1846— 0258 for two sections of data of ~2-years. Two cycles are shown for clarity. The solid line shows 
the pulse profile before the outburst in 2006 May. The dotted (red in the online version) line shows the summed pulse profile from 2008 
January to 2010 April. Subtracting the two profiles results in residuals with x^ u ~ 1-3, where the probability of this x^ v °r higher occurring 
by chance is 19%. Thus, the two profiles are not statistically significantly different. 
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Figure 4. Second derivative measurements for PSR J1846— 0258 from 2000 to 2010. The three measurements occurring after the glitch 
and outburst are one to two orders of magnitude larger than the subsequent measurements and vary in sign, indicating that they are 
severely contaminated by glitch recovery and/or timing noise. The inset shows measurements of ;> on a smaller scale to highlight the 
smaller variation away from the glitch recovery. The variation in u during the period from 2008 to 2010 is larger than in the pre-outburst 
era, and the mean value is systematically smaller. The one value of v pre-burst that is significantly larger than the average is directly after 
the candidate glitch near MJD 52910, and is the best evidence that a glitch actually occurred at that epoch. 
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Figure 5. Quantification of tlie timing noise in PSR J1846— 0258 over 10-year. Eacli point is a measurement of the A ■■■ parameter for 
approximately 2.5x 10^ s. This provides an estimate of the amount of timing noise observed in the pulsar, and shows a dramatic increase 
after the large glitch observed in 2006. 




Figure 6. Exposure-corrected Chandra ACIS images of the PWN associated with PSR J1846— 0258 in 1 - 7 keV. The 6" wide box along 
a position angle of 27° (north through east) is analyzed in detail for each epoch in Figure [71 
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Figure 7. Zoom-in of Figure [6] showing the northern clump of the PWN (boxed region in Figure O in the 1 — 7 keV range. The bottom 
right panel illustrates the count profiles extracted from the 6" wide box (as shown in Figure (6)1, indicating the evolution of the clump 
between epochs. 
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Table 1 

2008-2010 Spin Parmctcrs for PSR J1846-0258 



Phasc-cohcrcnt timing analysis 



Date range (Modified JuUan Day) 54492.089 - 55308.598 

Date range (Years) 2008 Jan 27 - 2010 Apr 22 

Number of TOAs 100 

Epoch (Year) 2009 Mar 1 

Epoch (MJD) 54834.0 

V (Hz) 3.0621185502(4) 

V (10-11 s-2) -6.664350(2) 

V (10-21 s"3) 2.725(3) 
Number of derivatives fitted 2 

RMS residuals (ms) 63^6 

Partially phasc-cohcrcnt timing analysis 

v (10-21 s-3) 3.13(19) 

Braking index in) 2.16(13) 

Quoted uncertainties arc the formal \u uncertainties as repored by TEMPO 
for the fully phase-coherent timing analysis. Details about the uncertainty on v and 
n for the partially coherent timing solution arc from a bootstrap analysis, as explained in the text (see Section 3.1). 



Table 2 

Variation of n with Number of Fitted Frequency Derivatives 



Derivatives 


Braking Index 


^i 


fitted 


n 




2 


1.888(2) 


39.71 


3 


2.010(3) 


6.12 


4 


1.980(8) 


6.00 


5 


2.05(1) 


5.48 


6 


2.08(2) 


5.50 


7 


2.51(4) 


3.51 


8 


2.60(4) 


3.26 


9 


2.95(7) 


2.93 


10 


2.91(7) 


2.73 


11 


2.1(1) 


2.12 


12 


2.1(1) 


2.14 



Braking index variation with number of fitted derivatives. The number 

of degrees of freedom for two fitted derivatives is 96. Uncertainties are the 

formal \a uncertainties returned by TEMPO for the timing solution 

spanning MJD 54492 - 55308. 



Table 3 

SPECTRAL FITS TO PSR J1846-0258 

T-i T aT t^ j-abs-FL /•unabs.Jr'ij j m /-abs-BB runabs.BtJ 2 / j r~ 

Epoch A^H r /o.s-io /0.5-16 ^^ /0.5-10 /0.5-16 X^/dof 

(10^2 cm-2) (10-12 erg (10-i2crg (keV) (10-^2 erg (10-i2crg 



2 s 1) cm 2 s 1) cm 2 s 1) cm 2 s i) 



2000 4^0^ 1.0+03 4.2 ±0.2 6.1 ±0.4 < 0.8 < 0.04 < 0.4 0.94/38 

2006 4.0t 1.9 ±0.1 1.3 ±0.3 31 ± 6 0.9 ± 0.2 1.7 ±0.2 3.2 ± 0.4 0.99/134 

2009 4.0t 1.0 ±0.1 3.6 ±0.5 5.2 ± 0.7 < 0.25 < 0.015 < 0.45 0.84/72 

t - held fixed in the fit. 

The uncertainties quoted are 90% confidence intervals and upper limits are at the 90% confidence level. 

The 2000 and 2006 results are from Ng et al. (2008). 



